Optical properties of random alloys : A formulation 
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We present here a formulation for the calculation of the configuration-averaged optical conductiv- 
ity in random alloys. Our formulation is based on the augmented-space theorem introduced by one 
of us [A. Mookerjee, J. Phys. C; Solid State Phys. 6, 1340 (1973)]. We show that disorder scattering 
renormalizes the electron and hole propagators as well as the transition amplitude. The corrections 
to the transition amplitude have been shown to be related to the self-energy of the propagators and 
vertex corrections. 

PACS numbers: 71.23.-k 



I. INTRODUCTION 

The object of our present study is to present a formu- 
lation for obtaining the configuration-averaged optical 
conductivity for random alloys. Because of randomness, 
there is a need to go beyond the usual reciprocal-space- 
based formulations for crystalline compounds. Instead of 
labelling the electronic states by the Bloch wave vector 
and band index (k, j), which is suitable for crystalline 
compounds, we have to label them by energy and the 
composite angular momentum L = {£,m,ms). In cases 
where the disorder is substitutional and homogeneous, 
in the sense that the occupation probabilities of lattice 
sites by atom species are independent of the site label, 
we can still label the configuration-averaged quantities by 
the reciprocal wave vectors. However, the band picture 
breaks down, and we cannot use the band index labeling 
of quantum states as in crystalline materials. Substi- 
tutional disorder dictates that we begin with a purely 
real-space representation and we have chosen the min- 
imal basis set of the tight-binding linearized muffin-tin 
orbitals (TB-LMTO) method [1, 2]. Configuration av- 
eraging over various random atomic arrangements has 
been carried out using the augmented-space formalism 
(ASF) introduced by us earher for the study of electronic 
properties of disordered systems [3, 7]. The ASF goes 
beyond the usual mean-field approaches and takes into 
account configuration fluctuations. This formalism has 
been described in detail in the referenced articles and 
the interested reader can go into the details in them. 
The contribution of the paper will be to show that 
the disorder-induced corrections to the averaged current 
terms in the optical conductivity are directly related ei- 
ther to the disorder-induced self-energy in the propaga- 
tors or to vertex corrections. Since the self-energy and 
the vertex corrections can be calculated for realistic bi- 
nary alloy systems, either within an augmented-space re- 
cursion [5-7] or within one of the mean-field approaches 



[4], this formulation will form the basis of subsequent 
calculations in real alloys. 



II. THE OPTICAL CONDUCTIVITY 

The Hamiltonian describing the effect of an external ra- 
diation field on the electronic states of a solid is given 
by 



^ ri 

1=1 



Here A(rj,t) and (j){ri,t) are the vector and scalar po- 
tentials seen by the zth electron because of the radiation 
field. We have used atomic units in which the electronic 
charge, mass, the Planck constant, and the velocity of 
light are set to unity {e—l, m—1, h=l, c=l). There are 
N electrons whose positions are labelled by . The effec- 
tive potential V(ri) experienced by individual electrons 
is expressed in the local-density approximation (LDA) of 
the density- functional theory (DFT). To first order in the 
vector potential, the Hamiltonian becomes 



^ n 1 

^ = E + ^('■') + A(r„i).j, . (1 

2=1 ' 



The current operator is related to the momentum oper- 
ator as (e/m) pj. In atomic units, = = Vj = dvi/dt. 
We work in the Coulomb gauge where V • A(ri,t) = 
and (j){ri,t) = 0, so that the electric field 



E(r„i) = - 



dt ■ 
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The Kubo formula relates the linear current response to 
the radiation field, 
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where t = t — t' and 6(t) is the Heaviside step function, 

. / 1 ifr > 
= |o ifr < 0. 

|$o) is the ground state of the unperturbed system, that 
is, the electrons in the sohd in the absence of the radia- 
tion field. In the absence of the radiation field, there is 
no photocurrent, i.e., (<f'ob/j|4'o) = 0. The fluctuation- 
dissipation theorem relates the imaginary part of the gen- 
eralized susceptibility to the correlation function as fol- 
lows: 



then the above equation becomes 

^(^) = li^E JdE [S^{E-,E++u,)... 

■■■+S^iE+,E- +aj)-S^{E+,E+ +uj) ... 
■■■ -S-,{E-,E- (5) 

where 

= \im f{E±iS). 

We have used the Herglotz property of the Green opera- 
tor. 



XU^) = ^(l-e-^-)5^.H, (2) 

where /? = l/fcsT; ks is the Boltzmann constant, T the 
temperature, and 



/oo 
dte"-Xf.Ar), z=a;+^0+ 
-OO 



and 



/>oo 
dt e'^^ ($o|Jm(t) >(0)|$o), z=w+^0+. 
-oo 

An expression for the correlation function, at T = K, 
can be obtained via the Kubo-Greenwood expression, 



dE Tr 



j, Im{G''(£)}jt 

ImiCiE + uj)}] . 



(3) 



We have assumed isotropy of the response so that the 
tensor 5^,^ is diagonal and we have defined <S'(ci;) as the 
direction averaged quantity i ^^^(ti;). j^y is • j, and 
e-y is the direction of polarization of the incoming photon. 
The operators G'"{E) and G''{E) are the resolvents of 
the Hamiltonian projected onto the subspaces spanned 
by the occupied and the unoccupied one-electron states. 

The trace is invariant in different representations. For 

crystalline systems, usually the Bloch basis {|k, j)} is 
used. For disordered systems, prior to configuration av- 
eraging, it is more convenient to use the real-space based 
screened (or tight-binding) mufHn-tin orbitals as a basis 

mm- 

If we define 



(4) 



G{E + i6) = G^'iE) - i sgn((5) G\E). 

For disordered materials, we shall be interested in obtain- 
ing the configuration averaged response functions. This 
will require the configuration averaging of quantities such 
as <S-y(2;i, Z2). 



III. CONFIGURATION AVERAGING 

Any description of a disordered system must be from a 
statistical point of view, since the properties of these sys- 
tems are random variables and any particular conGgura- 
tion is of little interest [38]. Consequently, the study of 
configuration-averaged properties of disordered systems 
has received much attention. 

Configuration averaging for response functions in disor- 
dered materials has had some history. The Ziman-Faber 
theory [8-10], much in use for liquids, is valid for elec- 
trons, weakly scattered from a dilute distribution of im- 
purities. The extended version of this theory was pro- 
posed by Evans et al. [11], but this too overlooks multiple 
scattering effects, as pointed out by Roth and Singh [12]. 
An effective medium approach (EMA) was proposed by 
Roth [13] and developed further by Roth and Singh [12- 
14] and Asano and Yonezawa [15]. This approach does 
take into account multiple scattering effectively and, as 
we shall sec, will have close similarities with the approach 
we propose in this paper. Velicky [16] has developed an 
expression for configuration-averaged response functions 
in random alloys for simple tight-binding models with 
one orbital per site and diagonal disordered within the 
coherent-potential approximation (CPA). Brouers and 
Vedyayev [17] have extended the formalism to transition- 
noble metal alloys. The CPA-like mean-field approach 
has been applied to response functions by Niizeki and 
co-workers [18-20] , who extended the pioneering work of 
Velicky [16] to longer-ranged random potentials. Mook- 
erjee [3] has introduced the ASF to tackle configuration 
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averaging. Within this formahsm he studied the role of 
macroscopic conservation laws on the response functions, 
leading to a Ward identity between the vertex corrections 
and the self-energy [21]. Within the CPA, vertex correc- 
tions were obtained by ingenious diagram summations by 
Leath [22]. There have been CPA calculations by Harris 
and Plischke [23] and Nauciel-Bloch and Riedinger [24]. 
In a series of papers, Mookerjee and co-workers [25-27] 
have applied the ASF to conductivity and optical conduc- 
tivity in random alloys. This will form the background 
of our present development. 

The ASF has been described thoroughly in a series of 
articles [3-7]. Here we shall introduce the salient fea- 
tures which will be required by us subsequently in this 
paper. We shall start from a first-principles tight-binding 
linearized muffin-tin orbitals (TB-LMTO) [1, 2] method. 

The quenched local randomness in the alloy is described 

by a set of random occupa tion variables ur, which takes 
the value 1 if the mufHn-tin labelled by {R} is occupied 
by an A- type atom and if it is occupied by a S. The 
atom sitting at {R} can either be of the type A {riR = 1) 
with probability xa or B [nu = 0) with probability xb- 

The ASF now introduces the configuration space $ of the 

set of N random variables {n/v} of rank 2^ spanned by 
configurations of the kind [ tti ••■ it •••), where 

I U) - \Ar) + \Br) 

I in) - \Aji) - \Br). 

If we define the configuration ] |t • • • | • • • ) as the average 
or reference configuration, then any other configuration 
may be uniquely labelled by the cardinality sequence, 
{Rk}, which is the sequence of positions where we have a 
J, configuration. The cardinaJity sequence of the reference 
configuration is the null sequence {0}. 



The augmented-space theorem [3] states that 

« A{{nn}) » = < {0}1A|{0} >, (6) 

where 



A({Mfi}) = J ■ J A{{Xr}) YldPiXR) e $. 

P{Xr) is the spectral density of the self-adjoint operator 
M.R, which is such that the probability density of hr is 
given by 

p{nR) = -- lim Im (t [ ({nR + i5)I - Mr) ^ \ t). 

Mij is an operator in the space of configurations iJjr of 
the variable ur. This is of rank 2 and is spanned by the 
states {[ Tk), I iii)}, 

Mr = XAV]i + XB V^R + V^Ixi" (r^J- + r^T) , (7) 

where = \ U){U \, K = I ii^){iR I and 

7^^ = I ti?)(ii? 1 are projection and transfer operators in 
configuration-space. 

Within the ASF, the configuration-averaged Green func- 
tion is given by [29] : 

« G(^) > = (1] (0I - H*^*^)"' ]1), (8) 



ff*^ = ^i Ap«o2: + BPflOP^ + FPfl(^(r^'^+r^T^ +^ Yrr, TRR^^ii, 

R [. R' ) 

(9) 

where 

A = Al6ll' = A{Cl/Al)/A{1/Al) 5ll', 

B = Bl6ll' = B[{Cl-z)/Al]/A{1/Al) Sll', 

F = FlSll' = F[{Cl-z)/Al]/A{1/Al) 6ll', 

V = Vll'{R-R') = A{1/Al)-^'^ Sll'{R-R') A{1IAl')-^'^ . (10) 



Vr = |-R)(-R| and Trri = \R){R'\ are projection and angular momentum index {f, m, mg}. C, A are the TB- 
transfer operators in real-space, and i is a composite LMTO potential parameters and S is the structure ma- 
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trix, in the most tight-binding (a) representation [1, 2]. 
The following functions are : 

A{y) = XAVA+XBUB, i.e., the average of y, 
B{y) = {xB - xA){yA- ys), 
F{y) = ^XBXA{yA -yB), 



and the contraction 

b{x'y b{x)^- = ie{t'-t) Srr^, 

where 

b{x) = bi{{t) = exp (ittlo^ buexp ^—itHo 
so that 



\1} = A m}+F |1) = 




We can reformulate the above in a second quantized for- 
malism. This follows the ideas put forward by Schultz 
and Shapero [28], which were extended in the ASF by 
Mookerjee [25]. For the real-space part, this is straight- 
forward, with a vacuum state described as that one which 
contains no electron-like excitations, and the fermion cre- 
ation and annihilation operators are a^j^^ and qhl for 
electrons at the site R with angular momentum indeces 
L. For the configuration-space part, we shall follow the 
ideas of Ref. [28] and consider the reference state to be the 
vacuum. Each spin Hip [39] at any site from up to down is 
then a creation of a configuration fluctuation. Since each 
site can have only two configurations, two up to down 
spin Hips cannot take place at a site. These excitations 
are then JocaJ and fermion-like. Each spin flip from down 
to up is a destruction of such a local pscudo-fcrmion. The 
Fock space is then spanned by all configuration states la- 
belled by the cardinality sequences. The corresponding 
fermion-like creation and annihilation operators are b\^ 
and These create and annihilate configuration fluc- 
tuations over the reference state. We should note that 
the configuration fluctuations are local and quenched. In 
second quantized form, the Hamiltonian becomes 

H = Ho -|- Hi, 

Ho = ^ al^^ttRL + X] ^RL,R'L' O^rl'^R'L', 

RL RL R'L' 

^1 = ^{^L a'RLO'RLb''RbR 
RL 

+ h a^RL^RL {bR + b^R)} , (fl) 

where 

{bR,bR'} = = {blb^], 



^{x,x') = ^im\{Tb{x')bHx)}m) 

= -ie{t' -t) 5RR,5{t-t'). (12) 

This pseudo-fermion formalism has been described earlier 
by Mookerjee [25] . The readers are referred to that article 
for further details. 



IV. AVERAGED GREEN FUNCTION IN THE 
PSEUDO-FERMION FORMALISM 

In this section, we shall develop a multiple scattering for- 
malism for the configuration-averaged Green function for 
a random binary alloy. The scattering is by configuration 
fluctuations and within the second-quantized formalism 
just described, the scattering diagrams are Feynman di- 
agrams. The formalism is very close to the Yonezawa- 
Matsubara scattering diagrams [15] and one can establish 
a one-to-one correspondence between them in the special 
case of diagonal disorder. 

The augmented-space theorem then states that 
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FIG. 1: The scattering diagrams for (a) the two topologically 
identical diagrams for n = 2 and (b) one of the 3! topologi- 
cally identical diagrams for n = 3. 
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< G{x,x') » = -i ^ 

n=0 



dtidt2 ■ ■ ■ dtn- 



{THi(ii)...H„(i„) a{x)a'<{x')}\0) 




FIG. 2: The topologically distinct scattering diagrams for 
n = 4. Top and middle are non-separable, the bottom is 
separable. The middle diagram is a skeleton one. 



where it is understood that the boldface operators are 
expressed by the matrix representation in {L} space and 

a{x) = aRL{t) = exp (^«tHo^ ttiiz, exp ^-iiHo^ , 
S = U(oo,-oo) and |0) = |O(g){0}), 

and 

u(i,i') = I- r dt" Hi{t")ij{t",t'). 

J — CO 

We may now apply Wick's theorem and Feynman's rules 
and generate a diagrammatic expansion for the averaged 
Green fimction <C G(a;, x') ^ in terms of the VGA Green 
function, 

g{x,x') = -I (0|{Ta(x)at(x')}|0). 
Let us examine a few terms in the series, 
(i) For n = 1, the term in the expansion is 

-i j dh (0| {THi(ti)a(a;)at(a;')} |0)conn. 

The contribution of this term is zero, since all 
three terms arising out of Hi [see eq. (11)] vanish 

because ({0}|6t(x)6(a;)|{0}), {{<l>}\b^x)\{<D}) , and 
({0}|6(a;)|{0}) are ah zero. 
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FIG. 3: Skeleton diagrams for order n = 5. 
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FIG. 4: The skeleton diagrams for the self energy S(z). 
(ii) For n = 2, the only nonvanishing terms come from 

J dh J dt2(0|{THi(ii)Hi(i2)a(x)a+(a;')}|0)conn 

= y dxi j dx2 (0| {Ta\xr)a{xi) . . . 
. . . a'^ {x2)a{x2)a{x)a^x')} |0)conn • • • 

. . . ({0}|r (&t(a;i) + (6t(x2) + b{x2)) |{0})con 

g{x,xi) g(xi,a:;2) g{x2,x') 7(.Ti,.-r2) ... 
+ s{x,X2) g{x2,xi) g{xi,x') 7(a;2,a;i)]. 



Figure 1 (a) shows a pictorial representation of the 
two terms, which are topologically identical and 
therefore have identical contributions. This cancels 
the (1/2!) term in the expansion for <C G(x, x') 
The F vertex [see eq. (11)] has a contribution F^i^r 
which is diagonal in L space Fl^/ = F^ 5ll', 
where 
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Cl 
Af 







(a^ 


Af). 



(iii) Figure 1 (b) shows one of the topologically identical 

diagrams (there are 3! = 6 such diagrams) for n = 
3. Note that it involves the scattering vertex B. 
This arises from the first term in the expression for 
Hi in eq. (11). Its contribution is also diagonal in 
L space Bl^i = B^ Sll', where 



_ {xb - Xa) 

Bl — , ^L- 

\/XbXa 

This scattering vertex cannot sit either in the left- 
most or in the rightmost positions, because one 
of the associated pseudo-fermion propagators van- 
ishes. 

(iv) For n = 4, there are two topologically distinct non- 
separable diagrams [40] : the double tent and the 
crossed tent diagrams shown as the two top dia- 
grams in Fig. 2. The inner tent in the top diagram 
goes on to rcnormalizc the interior Green function 
from g(a;, a;') to <C G(x,a;') !». As such, the mid- 
dle diagram is the only skeleton diagram at this 
order. There is one separable diagram (the bottom 
diagram in Fig. 2). This can be broken into two, as 
shown, by the dotted line. The rightmost diagram 
renormalizes the rightmost electron line. 

(v) The non-separable topologically distinct diagrams 
for n = 5 are shown in Fig. 3. We note that all 
odd-order diagrams must have an odd number of 
B vertices. 

The skeleton diagrams provide the expression for the self- 
energy for the Dyson equation which follows : 



< G(.x\.t') > = g{x,x')... 

■■■ + J dy jdy' g{x, y) S(2/, y') < G{y' , x') > . 

For homogeneous disorder, we have shown earlier that we 
have translational symmetry in the full augmentcd-space 
[29] . We can then take the Fourier transform of the above 
equation to get 



G{\<i,E) = g(k,i;) + g(k,i;) s(k,i;) G(k,i;). (13) 

The diagrams for the self-energy are shown in Fig. 4. In 
the above equation, each term is a matrix in {L} space. 
V. AVERAGED OPTICAL CONDUCTIVITY IN 
THE PSEUDO-FERMION FORMALISM 

We now go back to eq. (4) and discuss the configuration 

averaging of the two-particle Green functions of the kind 
S-^{zi,Z2)- The augmented-space theorem immediately 
implies that 



< 5^(zi,Z2) > 

= Tr ({0} fj^G"(^i)j;G^(z2)l {0} 



(14) 



The first thing to note about eq. (14) is that the right 

hand side is an average of four random functions whose 
fluctuations are correlated. The average of the prod- 
uct then involves the product of the averages and other 
contributions which come from averages taken in pairs, 
triplets, and all four random functions. 



Disorder-linduced renormalization of the 
current terms 



At this stage, in order to simplify notation, we shall omit 
the L index. However, we have to remember that all 
terms labelled by indeces ii or k are matrices in {L} 
space, so the order of multiplication of various terms in 
the expression has to be preserved. We shall also omit 
the 7 index of the current term indicating the required 
projection onto a direction. If required, we can put them 
back in the final expression. 

The real-space representation of the random current op- 
erator, can take the values j;^;^, , , j]^^, , or j^^, with 
probabilities x\ , xaxb, xbxa, and x^, respectively. We 
may rewrite jrr' as 



JRR' 



,AA 

3rr' nR UR' 

BA 



JhI' riR (1 - 
tir) nw (1 



nR') + ... 

- nfl)(l - uri). 



Following the same augmented-space procedure as for the 
single-particle Green functions, we get 
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(10) 



(15) 



(4) ^ 



(11) 



(5) (6) 




FIG. 5: The scattering vertices associated with the random current terms. 



J = 



R R' - 

... + {xb- xa)^/xaxb a)v a/?/ 

. . . +{XB - XAf j'kR' «fl o-R' b^R'bR- b^'nbR + xaXb a^j flij' {(^Ij + (^Ij. + 



(15) 



where 



JkK' ^ '^A yfifi' iRR') Xb {Jrr' — iRR'j , 

•(2) _ -AA I -BB _ -AB _ -BA 
JRR' ~ JRR' + JRR' JRR' JRR'- 

The first term in Fig. 5 is the averaged current. The fig- 
ure shows the 15 different scattering vertices arising from 
terms in eq. (15). The rule for obtaining the diagrams 
for the correlation function Sry{zi, Z2) is as follows : Take 
any two current diagrams from Fig. 5 and two propaga- 

I 



tors and join them end to end. Now join the configuration 
fluctuation lines (shown as dashed arrows) in all possible 
ways. 

The dominant contribution comes from the diagram 
shown in Fig. 6. Here the two current terms are the 
averaged current, and all configuration- fluctuation dec- 
orations renormalize only the two electron propagators. 
In this diagram the bold propagators are the fully renor- 
malized electron propagators and the contribution of this 
term is 



/ ^ «j(k)><G''(k,^i)>«j(k)>t<G^(k,02)>. (16) 

JbZ fSTT 

I 

We now focus on the main correction terms to the ex- pression in eq. (16). These are the correction terms to 
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FIG. 6: The diagram for [< j >< X j X >] 
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Renormalized Currents 
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(4) 



-A) 



(5) 



(6) 



(7) 



(8) 



Self-energy 



a 

F G 



G F 



FIG. 7: Relation between the renormalized currents and the 
self-energy. 



the averaged current which are closely related to the self- 
energies. They arise from a set of diagrams in which no 
disorder propagator (shown as dashed lines) joins cither 
two electron propagators or two of the current lines di- 
rectly. These diagrams are made out of a left renormal- 
ized current diagram chosen out of the diagrams (l)-(4) 
in Fig. 7 and one right renormalized current diagram from 
(5)-(8) connected by two renormalized propagators : the 
bottom one being a valence and the top a conduction 
electron propagator. 

Let us now obtain expressions for the renormalized cur- 
rents. A careful look at the sclf-encrgy diagrams (see the 
bottom of Fig. 7 and the example diagram shown there) 
shows that all self-energy diagrams have the structure 

S(k,2) = F{z) *(k,2) F{z), 

where $(k, z) is the Fourier transform of 

^RR'iz) = GrrM Pi5i.(z) Gr,r,{z). 



While the contribution of the diagram labelled (1) in 
Fig. 7 is 



j(i)(k) F(zi) *(k,^i) F(^i), 



(17) 



where 



Fll'{z) = ^/xaxb 



fUz) 

Fll' = y/XAXB Sll', 



SlL', 



the expression for (17) becomes 

f\k) F(zi) F-i(zi) S(k,zi) = j«(k) f(zi) S(k,zi). 

The contributions of the other diagrams in the left col- 
umn of Fig. 7 are 

E(k,Z2)F-i(z2)F(z2)j(^'(k) = S(k,Z2) f(^2) j^^Hk), 
F(^i) j^'Hk) F-i(^i) S(k,zi) = jW(k) f(^i) S(k,zi), 
S(k,^2) F-i(22) j('Hk) F(z2) = S(k,^2) f(22) j(''(k). 

Prom the forms of F(^;) and F{z), we note that : 

j«(k) F(z) F(z) =F(z) j(i)(k) F{z) =j«,(k)/^,(z). 

Similarly, the contributions of the diagrams in the right 
column in Fig. 7 are 

S(k, Zi) f(zi) j(l)(k), j(l)(k) f(z2) S(k, Z2), 
S(k, Zl) f (Zl) (k), j(^) (k) f (Z2) S(k, Z2). 

Closely related to the above diagrams is a group of di- 
agrams which describe joint fluctuations of one current 
and two propagators. Two such diagrams labelled (9) 
and (10) in Fig. 8 can also be expressed in terms of the 
self-energy : 

S(k,Z2) f(^2) j^'Hk) i{zi) S(k,zi), 
S(k,zi) f(2i) j(2)(k) f{z2) S(k,Z2). 

If we now gather all the contributions from these dia- 
grams, we may define a renormalized current term as 
follows : 
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FIG. 8: Few more renormalized currents. 



J««f(k,0i,^2) =«j(k)>+2 



S(k,02) f(^2) j(^Hk) + j^'Hk) f{zi) S(k,^i) 

... + E(k,^2) f(^2)j('Hk) f(-2l) S(k,0i). 

I 

The contribution of these disorder renormahzed currents and propagators to the correlation function is 
I 



(18) 



«; 5(1) (2:1, 02) > = 



BZ 



87r3 



Tr 



P*(k,2l,^2) « ^(k, zi) > r"(k,2i,^2)^ « G'=(k,22) > 



eff/ 



(19) 



(1) 



(2) 



r 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



FIG. 9: The scattering diagrams associated with joint fluctu- 
ations of the random current terms and one propagator. 



Wc now turn to terms which involve joint fluctuations 
between the two current terms and one propagator. We 
shall show that the corrections due to these terms are 
also related to the self-energy. Such diagrams arc shown 
in Fig. 9. Contribution of these diagrams is given by 



Jb 



Tr 



j(i)(k) f{zi) S(k,0i) f(^i) jW(k)' < G(k,02) > . 



BZ 

.. + j«(k)^ {{Z2) S(k,^2) f(^2) j'^Hk) « G{k,z,) » 
I 



(20) 



These terms have a slightly different structure than those self-energy. 

shown in Fig. 7. However, they still depend only on the Intuitively, we expect these to be the dominant disorder 
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scattering correction to the averaged current. It is im- 
portant to note that this correction can be obtained from 
the self-energy and is therefore eminently computation- 
ally feasible in the case of realistic alloys, once we have a 
feasible method for obtaining the self-energy. 

There are other scattering diagrams which are not related 

to the self-energy, but rather to the vertex corrections. 
In these diagrams, a disorder line connects both the elec- 
tron propagators directly. We expect these corrections 
to be less dominant. For the sake of completeness, we 
shall indicate how to obtain them in the Appendix. We 
should note that since these corrections are related to 
the vertex corrections and we shall indicate how to ob- 
tain them within a ladder approximation, we need not 
sacrifice these terms in a calculation for a realistic alloy 
if we do not wish to do so. However, in most cases we 
expect their contribution to be relatively small. 



B. The vertex correction 



both the propagators directly. These lead to vertex cor- 
rections due to electron-electron and electron-hole corre- 
lated propagation. Figure 10 show a few of these dia- 
grams. In general, we obtain a Bethe-Salpeter equation 
for the averaged two-particle propagator. We shall con- 
sider only one special class of vertex diagrams in this 
paper, namely the scattering diagrams which are built 
out of repeated vertices shown on the first line of Fig. 10. 
These arc called the ladder diagrams and can be summed 
up to all orders. This is the disorder scattering version 
of the random-phase approximation (RPA) for electron- 
electron scattering. There is another form of diagrams 
shown on the second line of Fig. 10 with ladder inser- 
tions between the crossed vertices. These arc known as 
maximally crossed diagrams. These diagrams lead to the 
localization effect. 

Here we shall sum the ladder diagrams to all orders. The 

contribution of the ladder diagram shown as the second 
diagram on the top line in Fig. 11 is 



Wc shall now examine the scattering diagrams we have 
left out, namely those in which disorder lines connect 



X! X! X! X! X! X! '^'R5Le,RlL^ Gr^Li,R2L2{zi) GR2L2,R3L3izi) ■■ ■ 

R1R2 R5 L1L2 L3L4 L^Lq 

■ ■ ■ •^'R3L3,RiL4, GRiL4.,R2L5 (^2) GR^L5,R5Le (^2) 



where Grlr'L'{z) = « GR^R'L'iz) > and 



Wt=FL{Z2) 



RL" ,RL' 



... + B^"{Z2) GRL",RL'iZ2)\\ Fl'{zi). 

Homogeneity in augmented-space means that this is in- 
dependent of R and it allows us to take the Fourier trans- 
forms leading to 



/BZ 

. . W 



j^^ ^ G(k,^2)pff(k,^i,^2) G(k,^i) ... 

y^^l^ G(k',^l) J^^^{k',ZuZ2) G{k',Z2) 

= r{zuZ2) W f{zuZ2) (21) 



We define 

87r3 



G(k,02) J'^"(k,^i,02) G{k,zi)=T{zuZ2) 



BZ 



/ 

Jbz 



87r3 



G(k',^i) r«(k',^i,^2)^G(k',^2) = r(^i,^2). 



r 



Let us now look at the contribution of the set of ladder 
diagrams. Each one of them has the same structure as 
eq. (21). We may then sum up the series as follows. 

Let us define 



^LiL2 

1^3 



lLlizi,Z2) = J^^ Gl3lA^,Zi) GL2Ld^,Z2) 



^lIlI = W^3 SL3L, 



These super-matrices in {L} space are written as A and 
w. The full ladder vertex may now be written as 



A{Z1,Z2) = ^ -I- ^ A('2lj -22) ^ + ^ A(-2lj -22) ^ A(-2l, -22) ^+ • • • 

= ^(I - Uzi,Z2) ^)"^ (22) 



The ladder diagram vertex correction now can be written 
as 
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Scattering diagrams 



p- 



-d 



' ' ~ 
-6 d — ^-6- 



Vertex Corrections 



q P 

d b 



Ladder 



iUlaximally crossed 



Others 




FIG. 10: The scattering diagrams leading to vertex corrections. 



«5ladder(^l,^2)> = Tr ^ ^ T^J (^i , ^2) A^J;[;^ T^^ (^1, ^2) = TrT{zi,Z2)®T{zuZ2) K{ZI,Z2). 



(23) 



L1L2 L3L4 



r 



VI. COMMENTS AND CONCLUSION 

Starting from the pseudo-fermion picture in the aug- 
mented space method, we have obtained an expression 
for the configuration-averaged optical conductivity. The 
disorder scattering renormalizes both the electron propa- 
gators as well as the current terms. We have shown that 
the dominant corrections to the averaged current can be 
related to the self-energy. For the sake of completeness, 
we have also shown that the remaining correction terms 
are related to the vertex corrections. We have also in- 
dicated how to obtain the vertex corrections within the 
ladder approximation. Once we set up a computationally 
feasible technique for the computation of the self-energy 
and the ladder approximation to the vertex correction, all 
the correction terms can be easily obtained. In an earlier 
communication [37], we have suggested the augmented- 
space recursion as a feasible technique for obtaining S(k) 
and have applied it for obtaining the complex band struc- 
ture and density of states of a series of realistic metallic 
alloys, namely AgPd and AuFe and most recently NiPt 



among others. We propose to use that technique and the 
results derived here to obtain the configuration-averaged 
optical conductivity in disordered metallic alloys. We 
intend to study, through numerical calculations, the rel- 
ative importance of the contribution of the difi'erent cor- 
rection terms. 



APPENDIX A: CORRECTIONS TO THE 
CURRENT TERM RELATED TO THE VERTEX 
CORRECTIONS 



For the sake of completeness, we shall also indicate the 

contribution of those scattering diagrams to the current 
which cannot be directly related to the self-energy but 
rather to vertex corrections. 

These remaining diagrams arc shown in the left column 
of Fig. 12. These diagrams cannot be related to the self- 
energies, but rather to specific vertex correction diagrams 
between the two propagators. There are three categories 
of diagrams : ones that involve F vertices [labelled (11)- 
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Real Space Representation 
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FIG. 11: The ladder scattering diagrams for the vertex correction in real-space and reciprocal-space representations. 






— ^ 














^ 



(17) 



(18) 



FIG. 12: Renormalized currents (left column) derived from vertex corrections (right column). 



(13)], ones that involve B [labelled (14)-(16)] type ver- 
tices and those that involve both [labelled (17)-(18)]. 

The first category of diagrams (11)-(13) contributes the 
following : 



F(^i) F-i(^i) j(2)(k) F(^2) F-\z2) A(^)(0,k; z,.Z2) 

^ diagram (11), 
F{z,) F-i(^i) j(2)(k) F{Z2) F-'iz,) A(^)(k,0;zi,Z2) 

=^ diagram (12), 

F(zi) F-\z,) j(^)(k) F(Z2) F-\z2) A(^)(k,k;zi,Z2) 

=> diagram (13). 
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Inserting the expressions for F{z) and F, we get a total 
contribution, 

Jl = f(^i) j(2)(k) f(22)|A(^)(0,k,2i,^2) ... 

. . . + A(^)(k, 0, 01, ^2) + A'^^k, k, 02, ^1)} (Al) 

Here, the vertex correction term A^^^ involves only F- 
like vertices in all four legs. Similarly, for the other two 
sets of diagrams we get 

B(0i) B-\z,) j(i)(k) F-\z2) A(^)(k,k;0i,02) 

^ diagram (14), 

F-i(0i) j(2)(k) B{Z2) B-\z2) A(^)(k,k;0i,Z2) 

=> diagram (15), 

B(zi) B-i(zi) j(2)(k) B{Z2) B-\z2) A(^)(k,k;0i,02) 

=> diagram (16). 

The total contribution wiU be 

...f(0i) j(2)(k) f(02)}A(^)(k,k;0i,02). (A2) 
I 



« 5(3) (01, 02) > = 




Finally, Fig. 13 shows the diagrams with joint fluctua- 
tions of two current terms and two propagators. These 
are also built out of vertex corrections. Note that each of 
the six diagrams can be broken up into a left and right 
part. For the diagrams shown in Fig. 13 all the right 

I 



A^-^) involves only i3-like vertices in its left-hand side 
legs. Finally for the last two diagrams, 



B{zi) B-\zr) j(2)(k) F{Z2) F-Hz2) A(^^)(k,k;0i,02) 

diagram (17). 

F(0i) F-i(0i) j(2)(k) B{Z2) B-\z2) A(^^)(k,k;zi,Z2) 

diagram (18). 



Their contribution is 



J3 = 2 f(0i) j(3)(k) f{z2) A(^^)(k,k;0i,02). (A3) 
Collecting together terms 



AJ = Jl + J2 + J3. 



The contribution of these disorder renormalized currents 
and propagators to the correlation function is 



(A4) 

I 

parts arc the same. Thirty other similar diagrams can 
be produced by replacing the right part with the five dif- 
ferent left parts mirror-imaged. The contribution of these 
diagrams is then, if 



AJ(k, zi, Z2) « G''(k, 01) > AJ(k, 01, 02)^ « G'(k, 02) > . 



K(k,0i,02) = f(02) j('^(k) f(0i)+b(02) j(2)(k) f(0i)+f(02) j(2)(k) b(0i) + b(02) j(2)(k) b(0i). 



where 
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BZ 



87r3 



Tr K(k,2;i,2;2) (g) K^{k,Zi,Z2) A(k, k, 2:1, 02). 



(A5) 




FIG. 13: Some of the scattering diagrams associated with 
joint fluctuations of the random current terms and two prop- 
agators. 
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